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Abstract. 

We present the analysis of the radial gradients of stellar absorption lines in a sample of 32 bulges of edge-on spiral galaxies, 
spanning nearly the full Hubble sequence (from SO to Sc types), and a large range of velocity dispersion (from about 60 to 
300 km s"'). Different diagnostics such as index-index, gradient-gradient diagrams, and simple stellar population models are 
used to tackle the origin of the variation of the bulge stellar population. We find that the vast majority of bulges show older 
age, lower metallicity and higher [ff/Fe] in their outer regions than in their central parts. The radial gradients in [Fe/H] are 2 to 
3 times larger than in Log(age). The relation between gradient and bulge velocity dispersion is interpreted as a gradual build 
up of the gradient mean values and their dispersions from high to low velocity dispersion, rather than a pure correlation. The 
bulge effective radii and the Hubble type of the parent galaxies seem to play a more minor role in causing the observed spatial 
distributions. At a given velocity dispersion, bulges and ellipticals share common properties. 

Key words, galaxies: abundances, bulges, formation, evolution, kinematics and dynamics, stellar content 



1. Introduction 

Combining results of several recent studies of galaxies of dif- 
ferent Hubble types, it seems clear that a proper understanding 
of bulges of spiral galaxies is crucial to understanding galaxy 
formation in general. It has been realized that early-type (cD 
to E/SO) galaxies harbor a variety of luminosity profiles, which 
has been interpreted as being du e to a varying contribution of 
a disk component (Saglia et al. Il997t de Jong et al. .2004 ). 



Bulges have been envisaged as former ellipticals which re- 
acted to the presence of disk accreted later on (e.g., Barnes & 
White [l984) . Other studies have advocated for a "secular evo- 
lution" scenario to build bulges, gradually inflating them using 
disc material (e.g., Sheth et al. 2005 and references therein). 
For all such studies, one faces the fact that the light distri- 
bution of (large) galaxies is dominated by two components, 
bulge and disk (Allen et al. 2006h . While total disk magni- 
tudes vary by ~ 2 mag along the Hubble sequence (taking the 
faintest or brightest extremes of the distribution at each Hubble 
type), the bulge's ones vary by twice this amount (Simien & de 



* Based on observations collected at the European Southern 
Observatory, proposals number 58.A-0I92(A), 59.A-0774(A) and 
6I.A-0326(A) 



Vaucouleurs 1 1 986t de Jong Il996h . i.e., properties of bulges are 
key to our understanding of the nature of the Hubble sequence. 

With the exception of the Milky Way and M3 1 , in which we 
can resolve individual stars (Sarajedini & Jablonka,2005; Olsen 
et al. 120061) . studies of bulges have to deal with integrated prop- 
erties of stellar populations. Despite the prospect of yielding 
crucial information on galaxy formation and assembly history 
scenarios from analyses of bulge properties, bulges have re- 
ceived significantly less attention than E and SO galaxies. This 
is likely a direct consequence of the considerable challenge 
of avoiding possible disk light contamination when studying 
bulges. 

Spectroscopic studies of the metallicity and age o f the cen- 
tral parts of bulges were pioneered by Bica (119881) . During 
the following several years, clear evidence was provided for 
a central metallicity-luminosity (Z - L) relation for bulges. 
Several studies underlined the similarities between ellipticals 
and bulges (Jablonka, Martin & Arimoto 1996; Idiart, de 
Freitas Pacheco & Costa 1996), both with respect to the slope 
of the Z-L relation and to the similar apparent a-element over- 
abundances (relative to solar abundance ratios) in the two types 
of systems. The Hubble type of the parent galaxy appeared to 
have little influence on the stellar population properties of the 
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bulge, which instead scaled mainly with central velocity disper- 
sion and bulge luminosity (Jablonka, Martin & Arimoto Il996h . 

A few years later, possible distinctions were suggested be- 
tween bulge properties of late-type spirals vs. t hose o f early- 
type spiral s (Falc on-Barroso, Peletier & Balcells l2002l: Proctor 
& Sansom liooi) . As pointed out by Thomas & Davies ('2006') 
however, these apparent differences seem to disappear when 
the appropriate range of central velocity dispersion is consid- 
ered (in particular at low values). It is also likely true that the 
homogeneity of the samples and their sizes have an impact on 
the drawn conclusions. 

An important limitation to our understanding of bulge prop- 
erties is that virtually all previous spectroscopic studies of 
bulges only sampled the light in the central regions which can 
easily be significantly affected by disk light, especially for late- 
type spirals. 

Substantial progress would be enabled with spatially re- 
solved spectroscopy of bulges, so that radial gradients of their 
stellar population can be measured. Unfortunately, investiga- 
tions of such radial gradients using large surveys have so far 
only addressed early-type galaxies i.e., elliptical and lenticular 
galaxies, with only very modest and rare excursions in to the 
case of later type galaxie s (e.g ., Sansom, Pe ace & Dodd ll994l: 
Proctor, S ansom & Reid bOOOl Ganda et al. l2006t Moorthy & 
Holtzman l2006l) . 

In order to assess the extent of the similarities between el- 
lipticals and bulges of spiral galaxies and to evaluate to what 
extent disk material may be influencing the formation and evo- 
lution of bulges, we have undertaken a very deep spectroscopic 
survey of bulges in a large, homogeneous sample of spiral 
galaxies spanning most of the Hubble sequence. 

The galaxy sample, the data reduction procedures and 
the indices are p resented in detail in Gorgas, Jablonka & 



Goudfrooij (|2007[ hereafter Paper I). We present the analysis 



of the line index radial changes in the current paper. It allows 
one at last to establish (or discard) the existence of scaling re- 
lations and their intrinsic dispersion, and to shed some light 
on the origin of the spatial variation of the stellar populations 
within bulges. 

2. The sample 

We selected a sample of 32 genuine (or close to) edge-on spi- 
ral galaxies. As shown in Paper I, 26 galaxies are inclined 
by 90 degrees, 6 galaxies have lower inclinations (60, 67, 
75, 80 and 88 degrees). Galaxies in the northern hemisphere 
were selected from the Uppsala General Catalog (Il973h . while 
southern galaxies were selected from the ESO/Uppsala catalog 
(Lauberts ll982h . Hubble types range from SO to Sc, with the 
following frequencies: 7 SO, 4 SO/a, 2 Sa, 4 Sab, 9 Sb, 5 Sbc, 
and 1 Sc. Their radial velocities range from 550 to 6200 km/s. 
Due to the high inclination of the galaxies, a precise morpho- 
logical classification is difficult. Therefore, we assign an uncer- 
tainty of about one Hubble type, which is a fair representation 
of the catalog-to-catalog variations for a given galaxy. 

Part of the spectroscopic observations were conducted at 
the 2.5m Isaac Newton Telescope (INT) of the Isaac Newton 
Group of telescopes on the island of La Palma (Spain). For the 



southern galaxies we used the 3.6-m ESO telescope and 3.5- 
m New Technology Telescope (NTT), both at ESO La Silla 
observatory (Chile). The spectrograph slit was oriented along 

the minor axis of the bulges. 

The recent work of Moorthy & Holtzman ( 20061) is the 
closest to the present study. The size of their sample is sim- 
ilar (38 galaxies vs. 32) and they also took care to span a 
large range of Hubble types. However, Moorthy & Holtzman's 
sample is composed of galaxies with any inclination angle, 
and only about ten of them are really close to edge-on. This 
leads Moorthy & Holtzman to investigate the transition be- 
tween bulge and disk populations, whereas we totally focus on 
the intrinsic bulge properties. The depth of the observations is 
also different: our spectra reach bulge regions at galactocentric 
distances <: 3 times larger in radius than theirs. 

3. The line-strength indices 

Our analysis is based on the measurements of line-st rength 
indices in the Lick/I DS sy stem (e.g., Worthey et al. 1 19941 : 
Worthey & Ottaviani '1997*). For the sake of a robust deter- 
mination of the radial variations of the stellar population by 
overcoming the noise in some individual indices at faint sur- 
face brightness level (as in the bulge outskirts), we introduce a 
few new indices that are linear combinations of some classical 
Lick/IDS ones. They are listed below: 

<Fe> - Fe2 - (Fe5270-i-Fe5335)/2; the classical mean Fe in- 
dex. 

Fe3 = (Fe4383-HFe5270-HFe5335)/3; first introduced by 
Kuntschner (2000). 

Fe4 = (Fe4383-HFe5270-i-Fe5335+Fe5406)/4 ; a coi nposit ion 
of indices dominated by Fe fines (see e.g. Worthey 19981 for 
the chemical content of the individual indices). 
Fe5 = (Ca4455-hFe4531+Fe5015)/3; based on indices domi- 
nated by a mixture of metals, including a elements. 
Fe6 = (Ca4455-HFe4531-HFe5015+Fe5709+Fe5782)/3; an ex- 
tended "mixed" composite Fe index. It could only be measured 
in galaxies observed with the 3.6m ESO telescope (see below). 
Ha - (H5A+HyA)/2 and Hf - {H6f+Hyf)/2 ; two composite 
H indices. 

Balmer = (H5a+H7a+H/S)/3; a composite Balmer index, only 
useful for regions free of emission lines, since H/3 is included. 

In order to be consistent in our fitting procedures, we use 
the same units for all indices, converting the index / into the 
magnitude index /' by: 



/' = -2.5 log 



1 



1 " 

- V h/Mk 



where n is the number of indices that compose the total index / 
and Alk is the width (in A) of the central passband in which the 
index 4 is measured. Gradients in atomic indices are more lin- 
ear with galactocentric radius when working with /' than with 
/. Another advantage of this strategy is that the values of all the 
indices are comparable with one another, since we are normal- 
izing by the bandpass width. 

We end up with a grand total of 33 different indices. 
Because the observations were done on different telescopes and 
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Fig. 1. Radial variation of a sample of 20 spectral indices for NGC 3957. Both sides of the galaxy minor axis are represented by 
gray squares and black circles for positive and negative radii, respectively. The dust lane is identified with open symbols. Dubious 
data are indicated by a star. Indices have been measured at the resolution of the Lick/IDS system. 
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NGC 585 H6'^ 




CNl 


CN2 


Ca4227' 


G-band 






Fe4383' 


Ca4455' 


Fe4531' 


-0.0086 


0.0083 


-0.0216 


-0.0165 


-0.0304 


-0.0488 


0.0395 


0.0292 


-0.0208 


-0.0001 


-0.0176 


0.0245 


0.0255 


0.0197 


0.0198 


0.0340 


0.0136 


0.0135 


0.0216 


0.0170 


0.0167 


0.0132 


C4668' 


w 


Fe5015' 


Mgi 


Mg2 


Mgb' 


Fe5270' 


Fe5335' 


Fe5406' 


Fe5709' 


Fe5782' 


-0.0215 


-0.0459 


-0.0205 


-0.0180 


-0.0360 


-0.0189 


-0.0132 


-0.0238 


-0.0131 


-0.0101 


-0.0087 


0.0116 


0.0331 


0.0094 


0.0065 


0.0064 


0.0071 


0.0075 


0.0078 


0.0109 


0.0062 


0.0130 


Na5895' 


TiOl 


Ti02 


<Fe>' 


Fe3' 


Fe4' 


Fe5' 


Fe6' 






B aimer 


-0.1275 


-0.0013 


0.0056 


-0.0186 


-0.0202 


-0.0184 


-0.0133 


-0.0116 


0.0155 


0.0188 


-0.0341 


0.0194 


0.0045 


0.0053 


0.0061 


0.0078 


0.0067 


0.0072 


0.0053 


0.0162 


0.0195 


0.0308 



Table 1. For each of our sample galaxy, the 33 measured gradients (first fine) and their attached errors (second line) are presented. 
The full table is accessible along with the electronic version of this paper. 



spectrographs, the wavelength coverage of the spectra is not 
totally uniform among galaxies. In particular, indices redder 
than 5600 A (i.e., Fe5709, Fe5782, Na5895, TiOi, and TiOa) 
could be measured only for the 7 galaxies that were observed 
at the ESO 3.6m telescope. We refer the reader to Paper I for 
a complete description of the conversion of the indices to the 
Lick/IDS system. 

4. The Gradients 

4.1. Gradient Measurements 

For each galaxy, the final radial sampling (i.e., the number 
of bins along the minor-axis radius) is a function of signal- 
to-noise ratio (S/N). Spatial rows of the 2-D spectra were 
summed until a minimum signal-to-noise of 10 per spatial bin 
was reached. This procedure did not concern the central parts, 
where the S/N values were much higher even in a single spec- 
trum row. At the end of this procedure, a set of 1-D spectra 
was associated with each galaxy, along with their correspond- 
ing distances to the galaxy centre. 

Bulge minor axis effective radii (r^//) were derived for all 
bulges as explained in Paper I. As illustrated in Figure [T] and 
Figure|2] our spectra reach the bulge effective radius for all our 
sample galaxies. For 25% of the sample, we reach 3 r^ff and 
for another 25%, we could reach 10 r^ff. In order to render the 
radial gradients fully comparable from one galaxy to another, 
we scaled all galactocentric distances by these effective radii, 

/' =Alog(r/reff) + B 

. Radial gradients were derived by taking into account both 
sides of the galaxy minor axes, avoiding dubious measurements 
such as those placed at the edges of the CCD, excluding loca- 
tions in the dust lanes, locations affected by emission lines, as 
well as the inner 0.5 - 1 arcsec (depending on the slit width) 
as they could be aff'ected by disk light and seeing effects. The 
limits of the dust-obscured zone, translated in a decrease in lu- 
minosity along the slits, have delineated the edges of the dust 
lanes. 

We adopted an error-weighted least-squares linear fit to the 
data, with a two-iterations procedure to reject points falling be- 
yond 5cr from the mean relation. The gradients for our 33 in- 
dices and 32 galaxies is accessible with the electronic version 
of this paper. Table [Uprovides an example of its format. 



An example of the resulting gradients is presented in 
Figure [T] for NGC 3957 and a subset of 20 indices. The plain 
lines indicate the fitted gradients and the error bar at the left 
edge of the line shows the residual standard deviation of the fit. 
Figure|2]displays the radial profile of the Mg2 index for the full 
sample of galaxies and gives an idea of the variety of gradients. 
Galaxies like NGC 5084, which represents the most extreme 
case, flatten in the inner regions. For most galaxies, however, 
a radial linear variation of the indices with radius is a very fair 
representation. 

4.2. Gradient Amplitudes and Statistics 

Figure[3]presents the distribution of the gradient amplitudes for 
our main indices. The values of the atomic gradients are small 
relative to past studies due to the logarithmic scale we are us- 
ing. However, as seen in Figure [1] the indices can vary by 30 
to 50 % from the central to the outer bulge regions. The vast 
majority of the galaxies exhibit negative gradients in all indices 
(except for the ones involving Balmer lines), i.e., the absorption 
features increase in strength towards the centre of the bulges. 
The range of gradient amplitudes is quite similar from one in- 
dex to the other Another interesting feature of Figure [3]is that, 
despite a non-negligible dispersion, the distributions are nearly 
always well peaked, suggesting that some factor other than ran- 
domness is at the origin of the spatial distribution of the bulge 
stellar population properties. 

Table |2] gives the error-weighted mean gradients and the 
dispersion of their distributions for our 33 indices. A compar- 
ison between the mean values (column 2) and their attached 
errors (column 3) demonstrates that the negative values of the 
gradients are significant. The dispersions of gradients among 
galaxies are also significant and not due to observational errors 
as indicated by the intrinsic dispersions. This means that there 
is a variety of stellar populations in bulges and an intrinsic di- 
versity in their spatial variation, the origin of which we will try 
to unveil. 

The dispersion around the mean gradient amplitudes (ctj^js) 
decreases with redder indices (i.e., those measured in the red- 
der part of the bulge spectrum). While this is partly due to the 
measurement errors which get smaller at larger wavelength, 
this cannot fully explain the measured dispersions, as can be 
inferred from the intrinsic ones. The indices which measure a 
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Index 



Mean 



(2) (3) 



cr 
(4) 



(5) 



H6'^ 


0.0236 


0.0218 


0.0123 


0.0180 




0.0159 


0.0204 


0.0149 


0.0140 


CN, 


-0.0618 


0.0335 


0.0126 


0.0310 


CNz 


-0.0635 


0.0359 


0.0154 


0.0325 


Ca4227' 


-0.0084 


0.0148 


0.0156 




G4300' 


-0.0135 


0.0183 


0.0123 


0.0135 


HTa 


0.0339 


0.0174 


0.0095 


0.0146 




0.0314 


0.0226 


0.0111 


0.0197 


Fe4383' 


-0.0265 


0.0101 


0.0102 




Ca4455' 


-0.0162 


0.0125 


0.0106 


0.0067 


Fe4531' 


-0.0135 


0.0108 


0.0080 


0.0072 


C4668' 


-0.0295 


0.0170 


0.0072 


0.0154 


H/3' 


-0.0039 


0.0228 


0.0190 


0.0125 


Fe5015' 


-0.0158 


0.0083 


0.0064 


0.0053 


Mgi 


-0.0264 


0.0157 


0.0058 


0.0146 


Mg2 


-0.0459 


0.0205 


0.0058 


0.0197 


Mgb' 


-0.0223 


0.0127 


0.0071 


0.0105 


Fe5270' 


-0.0149 


0.0095 


0.0066 


0.0068 


Fe5335' 


-0.0166 


0.0087 


0.0067 


0.0055 


Fe5406' 


-0.0145 


0.0097 


0.0069 


0.0068 


Fe5709' 


-0.0070 


0.0064 


0.0080 




Fe5782' 


-0.0177 


0.0112 


0.0107 


0.0034 


Na5895' 


-0.0881 


0.0632 


0.0146 


0.0615 


TiOi 


-0.0025 


0.0094 


0.0060 


0.0072 


TiOj 


-0.0064 


0.0117 


0.0062 


0.0100 


<Fe>' 


-0.0163 


0.0069 


0.0051 


0.0047 


Fe3' 


-0.0198 


0.0075 


0.0052 


0.0054 


Fe4' 


-0.0186 


0.0070 


0.0043 


0.0056 


Fe5' 


-0.0146 


0.0076 


0.0052 


0.0056 


Fe6' 


-0.0099 


0.0052 


0.0069 






0.0293 


0.0177 


0.0077 


0.0160 




0.0239 


0.0212 


0.0093 


0.0191 


B aimer 


0.0202 


0.0180 


0.0157 


0.0088 



Table 2. Statistics of the gradients for our set of indices. 
Column (2): error-weighted mean values of gradients; (3): dis- 
persion around the mean; (4): expected dispersion from the er- 
rors; (5): intrinsic dispersion. 



given class of chemical elements (i.e., a elements, iron, CNO) 
present similar intrinsic dispersions. At this stage, it is difficult 
to identify the origin of the differences in intrinsic dispersion 
between the different indices. They can have two and possibly 
mixed origins: varying ratios of the different chemical species, 
but also, at least partly, difference in the dynamics of the in- 
dices, i.e. their relative sensitivity to changes in abundance or 
age. 



5. The Scaling relations 

As seen in the previous section, bulges, in general, do exhibit 
radial changes in their stellar population properties: their cen- 
tral indices are the strongest. In this section, we investigate 
whether the strength of the gradients scales with some struc- 
tural parameters or dynamical properties of the bulges. 



5.1. The existence of the gradient-o-Q relation 

Figure |4]presents the relation between the gradient amplitudes 
and the bulge central velocity dispersion (ctq, in km s"'), calcu- 
lated in a 2"-radius aperture. The solid lines indicate the error- 
weighted least-square linear fits: 

A/' 



A log(r/reff) 



■■ A log(cro) + B 



Table |3] provides the coefficients A and B, along with the 
error-weighted standard deviation cr. In order to evaluate the 
significance of the fits, we performed an F-test which checks 
the validity of the assumption of a description of the data by 
a flat relation fixed to the mean of the gradients. It compares 
the variance of the distribution implied by such an assumption 
to the observed one. The smaller the probability returned by 
the F-test, the more different the distributions. We also used 
the Spearman test, which qualifies the significance of the cor- 
relation between the gradients and ctq. This test does not take 
the observational errors into account. This is why we use two 
independent statistical tests. The smaller the Spearman proba- 
bility, the more significant the correlation between the gradi- 
ents and the bulge central velocity dispersions. As also shown 
in Figure m there seems to be no statistically significant vari- 
ation of the amplitude of the stellar population gradients with 
the bulge central velocity dispersion, with the exception of the 
Mgi and Mg2 indices. 

Intriguingly, while Mgj and Mg2 (both of which are sen- 
sitive to magnesium abundance) vary with erg, Mgb does not. 
This dichotomy between the behaviors of Mg2 and Mgb has al- 
ready been noticed for elliptical galaxies. Mehlert et al. (l2003b 
concentrate on Mgb gradients in elliptical galaxies with a range 
of velocity dispersions similar to ours and do not find any corre- 
lation. More recently, Sanchez-Blazquez et al. 12006) confirm 
this absence of correlation. Meanwhile, Carollo et al. (Il993h 
reported on the correlation of Mg2 gr adients with ctq among 
elliptical galaxies. Ogando et al. ( 2005b also observe a relation 
between the Mg2 gradients and the velocity dispersions of their 
galaxies. 

If one consi ders th e simple stellar p opulat ion (SSP) mod- 
els of Vazdekis (Il999h or Thomas et al. (l2003h . the sensitivity 
of Mgb' to age at fixed metallicity (or to metallicity at fixed 
age) is ~ half that of Mg2. Indeed, our Mgb' gradients are half 
as large as the Mg2 ones on average. This is however not suf- 
ficient to explain the change in behavior of the two indices. It 
has to be attributed to the intrinsic properties of the indices and, 
in particular, to their dependencies to abundances. As a clear 
illustration of this, four galaxies (UGC 10043, ESQ 443-042, 
ESQ 512-012, and IC 1970) do not exhibit any gradient in Mgi 
nor in Mg2, but they do have significantly negative Mgb' gra- 
dients. In those cases, C4668' shows no radial change either 
Therefore, for those galaxies, the carbon abundance seems to 
supersede the magnesium sensitivity of the indices Mg2 and 
Mgi. This Mg vs. C balance is seen in a quite spectacular way 
for these 4 galaxies, but it could play a (likely more subtle) role 
in other bulges. In any case, the sensitivity to different chemical 
elements seems to be able to explain the discrepant behavior of 
the three "magnesium" indices. 




Fig. 4. The relation between the gradients and the bulge central velocity dispersions for the 25 indices that are available for the 
whole sample of bulges. The color code for the different Hubble types is the following: Black corresponds to E/SO, red to SO/a, 
orange to Sa, magenta to Sab, green to Sb, blue to Sbc, and cyan to Sc types. The black solid lines show the linear fits. 
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Fig. 5. Left panels: The relation between Mg2 gradient amplitudes and central velocity dispersion for bulges in this paper (lower 
left panel) and for elliptical galaxies in the sample of Sanchez-Blazquez et al. (2006; upper left panel). To guide the eye, a 
least-squares fit to the data for the bulges in this paper is plotted as dotted lines in both left panels. Right panel: Histograms of 
Mg2 gradient amplitudes for bulges with central velocity dispersions cr < 125 km s"' (dashed lines and grey hashing) and with 
cr > 125 km s"' (solid lines). 



Index 





0.025+ 0.021 


-0.032+ 0.047 


0.02 


0.97 


0.81 




0.028± 0.025 


-0.045+ 0.055 


0.02 


0.95 


0.98 


CNl 


-0.046± 0.022 


0.039+ 0.049 


0.03 


0.95 


0.17 


CN2 


-0.033+ 0.026 


0.010+ 0.057 


0.04 


0.98 


0.11 


Ca4227' 


-0.058+ 0.026 


0.117+ 0.056 


0.02 


0.56 


0.20 


G4300' 


-0.036+ 0.020 


0.065+ 0.044 


0.02 


0.89 


0.02 




0.042+ 0.015 


-0.059+ 0.034 


0.02 


0.83 


0.10 


Wf 


0.042+ 0.017 


-0.061+ 0.038 


0.03 


0.89 


0.05 


Fe4383' 


-0.015+ 0.015 


0.008+ 0.034 


0.01 


0.92 


0.37 


Ca4455' 


-0.005+ 0.017 


-0.006+0.036 


0.01 


1.00 


0.65 


Fe4531' 


-0.021+ 0.013 


0.032+ 0.028 


0.01 


0.89 


0.06 


C4668' 


-0.009+ 0.012 


-0.009+0.026 


0.02 


0.99 


0.53 


Fe5015' 


-0.003+ 0.011 


-0.009+ 0.025 


0.01 


1.00 


0.08 


Mgl 


-0.066+ 0.010 


0.121+ 0.021 


0.01 


0.49 


0.01 


Mg2 


-0.053+ 0.009 


0.071+ 0.020 


0.02 


0.79 


0.00 


Mg; 


-0.004+ 0.012 


-0.014+ 0.026 


0.01 


1.00 


0.42 


Fe5270' 


-0.009+ 0.010 


0.006+ 0.023 


0.01 


0.97 


0.33 


Fe5335' 


-0.008+ 0.011 


0.002+ 0.024 


0.01 


0.97 


0.15 


Fe5406' 


-0.019+ 0.010 


0.027+ 0.023 


0.01 


0.87 


0.08 


<Fe>' 


-0.014+ 0.008 


0.015+ 0.018 


0.01 


0.88 


0.13 


Fe3' 


-0.014+ 0.008 


0.011+ 0.018 


0.01 


0.90 


0.15 


Fe4' 


-0.012+ 0.007 


0.008+ 0.014 


0.01 


0.90 


0.07 


Fe5' 


-0.007+ 0.008 


0.002+ 0.018 


0.01 


0.97 


0.06 




0.026+ 0.013 


-0.027+ 0.028 


0.02 


0.94 


0.07 




0.036+ 0.015 


-0.054+ 0.033 


0.02 


0.92 


0.16 



Table 3. Linear fits of the variation of the gradients with the 
bulge central velocity dispersion for the full galaxy sample. A 
stands for the slope of the relations, B for their zero points, and 
cr is the dispersion of the fits. F and S give the probabilities of 
the F- and Spearman , respectively (see text). 



The indices primarily sensitive to iron exhibit gradients 
whose amplitudes seem independent of the bulge velocity dis- 
persion. As can be seen in Table |2l the intrinsic dispersions 
of these amplitudes are extremely small and leave hardly any 
room for any variation beyond the one consistent with the un- 
certainties. 

5.2. Building- up the gradient-o-Q relation 

The galaxy Hubble types are color-coded in Figure HI All 
galaxy types span nearly the full range of central velocity dis- 
persions, with the exception of the earliest-type galaxies (E/SO 
and so/a) whose central velocity dispersions do not reach low 
values. 

For the E/SO to SO/a galaxies in our sample, the measure- 
ment errors are extremely small and therefore give an excellent 
idea of the intrinsic dispersion at fixed velocity dispersion. We 
find that their gradients are independent of the bulge velocity 
dispersion, whatever index is considered. The range of Mg2 
gradients in our sample is entirely consiste nt with that of the 
compilation of SO galaxies in Ogando et al. ( 12005b . 

If one only considers galaxies of types later than Sa, then 
the G-band and Fe5406 indices reveal some variation with ctq 
as well as Mgi and Mg2, and the slopes of the relations of the 
latter two indices with ctq are slightly stronger than when the 
full sample is considered. Therefore, it looks like the Mgi and 
Mg2 trends seen in the full sample arise mainly from the Sa 
through Sc galaxy types. 

However, rather than interpreting the trends between gradi- 
ent amplitudes and cry as direct correlations, we find that they 
are due to the combination of a gradual and differential pop- 
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ulation of the diagrams together with the dynamical range of 
the indices (i.e., their intrinsic capacity to reach high values 
by nature). Rather than a smooth change in gradient amplitude 
with velocity dispersion, one witnesses a three-step process: (i) 
At large cr, the dispersion among gradients is large but small 
gradients are relatively rare, (ii) At smaller cr, the dispersion 
remains large, but galaxies with very weak gradients appear in 
larger number, (iii) Finally, at cr < 125 km s"', one finds only 
bulges with either very weak or no radial variation of their stel- 
lar population. These successive regimes are, of course, best 
seen with indices which have a large dynamical range of such 
as Mgi and Mg2, but it is likely present in all indices. 

Only a few previous studies on gradients included (early- 
type) galaxies with central velocities below ~ 125 km s 



(Gorgas et al. 119971: Sanchez-Blazquez et al. l2006 h. Just as 
found here, their low-cr galaxies do have very weak or null 
gradients. Gorgas et al. Jl997h find evidence for their bulges 
to have shallower gradients than ellipticals. However, they 
compare systems with very different velocity dispersions. The 
study of Ogando et al. (l2005b notices a decreasing number 
of E and SO galaxies harboring steep Mg2 gradients with de- 
creasing velocity dispersion, which is similar to what we find. 
Finally, Moorthy & Holtzman (l2006t) mention that the bulges 
of their sample showing insignificant gradients in [MgFe] have 
small sizes. Histograms of the distribution of Mg2 gradients 
for bulges with cr < 125 km s"' vs. those with cr > 125 km s"' 
are shown in Figure |5] Formally, a Kolmogorov-Smirnov test 
indicates that the two distributions are different at the 99.95% 
probability level. We also show a com parison with the ellipti- 
cal galaxies of Sanchez-Blazquez et al. (|2006l) which illustrates 
that Mg2 gradients disappear at low velocity dispersion in these 
systems as well. 

5.3. Other scaling relations 

Mehlert et al. (l2003h report on a stronger correlation of ab- 
sorption line strength gradients with the galaxy velocity dis- 
persion profiles than with their central values. We indeed find 
a correlation, but to a level which is no more significant than 
with cro. This result holds both when considering the full sam- 
ple and when we restrict our sample to SO and SO/a galaxies. 
This is in agreement with the results of Sanchez-Blazquez et al. 
(l2006h Ibr elliptical galaxies. Figure|6]displays the relations be- 
tween the bulge gradients and their effective radii. Here again, 
Mgi and Mg2 are the only indices for which a dependence is 
clearly seen. However, the slopes of the relations are 2 - 3 times 
smaller than the ones with cry. This holds for all indices. This 
means that any dependence on the bulge mass is dominated by 
the effect of the the bulge central velocity dispersion. 

6. Index-index diagrams 

We now try to understand the nature of the radial gradients 
in terms of physical quantities such as age and/or chemical 
abundance variations of the stellar population. We have con- 
sidered the solar-scaled SSP mo dels of Vazd ekis (il999 ) and 
the SSP models of Thomas et al. (I2003ll2004 which allow a- 
enhancement variations. The use of these two series of models 



was meant to assess the robustness of our conclusions. As a 
matter of fact, we find that they indeed yield the same broad 
conclusions, and we consequently choose to only present the 
analysis conducted with the Thomas et al. models which in- 
clude non-solar [a/Fe] values, and thus offer a more detailed 
description. 

One shortcoming using SSPs is that when directly applied 
to galactic composite stellar populations, the absolute values of 
age and metallicity are intrinsically biased, due to the simpli- 
fied description as "single" stellar populations. They are meant 
to get close enough to the luminosity-weighted mean quanti- 
ties. Another problem is the existence of some degeneracies, 
in particular between [a/Fe] and age for a fair number of in- 
dices. For these reasons, we avoid assigning absolute ages or 
metallicity values to our galaxies. 

As mentioned earlier, galaxies can roughly be divided in 
two categories: Those with significant radial gradients and 
those without. In the following we address each group sep- 
arately with the aim of possibly identify differences between 
them. 



6.1. Galaxies with strong gradients 




<Fe>' 

Fig. 7. Index-index diagram for galaxies harboring clear gradi- 
ents. SSP Model grids are shown for ages running from 15Gyr 
to 3Gyr The location of [Fe/H] = -0.33, 0.00, 0.35 and 0.67 
are identified by dotted lines. At the same time, they indicate 
the direction of an age variation. The gray and black grids cor- 
respond to [afFe]-Q.3 and 0.0, respectively. Arrows connect 
the central indices, integrated in a 4" aperture, to the values 
of the indices at the effective radius of the bulge. The galaxy 
Hubble type color coding is the same as in Figure H] 



First, we consider galaxies for which Mg2 and <Fe>' gradi- 
ents are strictly negative, taking into account their errors (Icr). 
Mg2 is a composite index, simultaneously sensitive to metal- 
licity, a-elements, and age. Therefore, this criterion does not 
select any particular type of galaxy. The criterion on <Fe>' dis- 
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cards 3 bulges, whose gradients, although negative, have large 
errors. Including them would not help our understanding. 
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Fig. 8. Index-index diagram for galaxies harboring clear Mg2 
gradients. The two grids of models for [afFe]-Q.Q and 0.3 are 
shown by black and light gray lines, respectively. Ages run 
from 15 Gyr to 3 Gyr and go along the dashed lines (younger 
ages correspond to higher H^). One Gyr separates two adjacent 
age lines. Arrows join the central indices, integrated in a 4 arc- 
sec aperture, with the value of the indices at the bulge effective 
radius. The galaxy Hubble type color coding is the same as in 
Figure |4] 

While the Mgb' vs <Fe>' diagram is highly degenerate in 
age and [Fe/H], it does offer a good determination of [a/Fe]. 
In Figure |7] for each galaxy, we join with an arrow the central 
indices to the indices measured at the bulge effective radius. All 
galaxies lie within a region delimited by the two lines at [a/Fe] 
= and [a/Fe] =0.3 and the two iso-metallicity lines at [Fe/H] 
values of -0.33 and 0.67. The iso-age lines span the range 3- 
15 Gyr. As shown by the direction of the arrows, most of the 
galaxies show both a decrease in metallicities and an increase 
in [of/Fe] towards larger radii. The variation in [Fe/H] is of the 
order of 0.3 dex, while it is smaller for [a/Fe] 0.1 dex). A 
few bulges seem to harbor pure [a/Fe] or pure [Fe/H] gradients. 

Figure |8] presents the variation of Fe4' as a function of 
H^ for our galaxies, superimposed on the same grids of mod- 
els as mentioned above. Here again, the index gradients are 
clearly mostly due to metallicity. H^ gradients are small and 
positive. An interesting feature is that the longest arrows (i.e., 
the largest gradients) are located at the bottom of the diagram, 
implying a slightly older luminosity-weighted mean age. The 
relative shift is small, of the order of ~2Gyr, but the effect 
seems systematic. Taking into account the increase in [a/Fe] 
at larger radii revealed by the previous diagnostic diagram, the 
amount of possible variation due to age corresponds to a few 
Gyrs at most (often of the order of 1-3 Gyr), the central parts 
of the bulges being younger than their outer regions. There are 
three apparent exceptions to this general trend: UGC 11552 
and UGC 11587 seem to have stronger age gradients than the 
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Fig. 9. Gradient-gradient diagrams for galaxies harboring clear 
Mg2 gradients. The galaxy Hubble type color coding is the 
same as in Figure |4] The direction of variation of the indices 
with [Fe/H] only, or with age only, is shown by the plain and 
dashed lines, respectively. 



bulk of the galaxies, as they have strong negative H^ gradi- 
ents. The third galaxy, IC 5176, harbors a strong and posi- 
tive H^ gradient. The bulge of this galaxy exhibits gradients 
in C4668' and in Mgb' that are among the smallest in the sam- 
ple, while showing a strong radial increase of its high-order 
Balmer lines towards the outer regions of the bulge. Its outer 
parts therefore seem indeed younger than the central ones. As 
a matter of fact. Ha and [N ii] emission are present all the way 
from the central to the outer regions, which increases the un- 
certainties of both Balm er and Mgb index measurements (e.g. 
Goudfrooij & Emselleml lT996!). Admittedly, we do not see any 
straightforward explanation for these traces of gas. One possi- 
bility is pollution by the disk co mponent, even though I C 5176 
does not (currently) have a bar dChung & Bureaull2004 . 
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Another way to look at the question of the nature of gradi- 
ents is provided by Figure|9l The SPP models, plotted in index- 
index diagrams, show grids from which one can easily derive 
the amount of variation when only one of the parameters (age, 
metallicity, or [a/Fe]) is varying. 

The advantage of this method is that it is immediately ap- 
plicable to any radial gradient-gradient observational plot. The 
slopes associated with expected pure age and pure metallicity 
variations have been derived for stellar populations of ages be- 
tween 3 and 15 Gyr and metallicities ([Fe/H]) between -1.35 
and 0.67. We considered [a/Fe] = 0.3, given our diagnosis of 
Figure |2l The results would change in detail but not in con- 
clusion if we had used [a/Fe] = 0.0. The slope of the metal- 
licity variation depends on the age considered and, conversely, 
the age variation depends on the range of metallicity consid- 
ered. For this reason, we show the maximum, mean and min- 
imum slopes. Figure |9] as was Figure [8] is meant to summa- 
rize our investigations. Other sets of indices have been consid- 
ered, however, we avoid any redundant illustration. The choice 
of the two sets of indices, (Mg2, C4668') and (H^, Fe4') was 
driven by four criteria: (i) They involve well-measured indices 
(ii) the index-index model grids are regular, (iii) they show sig- 
nificantly different slopes for age and metallicity changes, and 
(iv) the diagrams include sensitivities to all three parameters 
(age, [Fe/H] and [a/Fe]), so that the three can be considered 
together. 

Figure |9] demonstrates clearly that a pure metalUcity gra- 
dient nearly fully explains the amplitude of the observed gra- 
dients. Interestingly, the smallest slopes shown in Figure |9] fit 
the data the best (they correspond to the oldest ages). Only one 
galaxy, IC 5 176, is more compatible with a major radial change 
in age, as discussed above. Here the line passing through the 
IC 5176 location coiTesponds to metallicities below solar. 

All galaxies do not fall exactly on the models tracing a pure 
metallicity variation. The residual distances to these lines are 
compatible with an age and/or [a/Fe] variation. As these two 
parameters make the indices move in the same direction, it is 
very hard to distinguish between them using this method, but 
we already know that likely both can be accounted for, in small 
proportions. 

6.2. Galaxies with weak or no gradients 

We define this category of galaxies as bulges whose gradients 
in Mg2 are compatible with a null value within the errors. Eight 
galaxies in our sample fall in this category, i.e., 25 percent of 
our sample. Their Hubble types are mixed, but no very early- 
type galaxy is present (2 Sab, 2 Sb, 3 Sbc, 1 Sc). The three Sbc 
galaxies (UGC 10043, ESO 512-012, NGC 1351A) have large 
error bars, both in Mg2 gradients and in log cr, essentially due 
to the restricted number of bins on which the gradients could 
be evaluated robustly. 

As we did for the galaxies with strong gradients, we ex- 
amine the relation between Fe4' and H^ in Figure [TOl All the 
bulges with weak gradients are shifted to lower metallicities, as 
compared to the strong gradient ones. 
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Fig. 10. Index-index diagram for galaxies harboring weak or no 
Mg2 gradient. Ages run from 1 Gyr to 15 Gyr. The rest of the 
symbols and parameters are the same as in Figures[8]and|4] 

Furthermore, none of these bulges show strong Fe4' and 
H^. This means that weak gradients are associated with low 
global chemical enrichment. More precisely, the indices mea- 
sured at Teff populate a region of the diagram which is very 
similar to that of the bulges with strong gradients (cf. FigureO, 
but the central index strengths do differ significantly. This sup- 
ports a scenario where star formation proceeds from the outer 
to the inner parts, and where the outer bulge regions share 
universal chemical properties. The diff'erence between bulges 
arises from star formation in their central regions. 

Two galaxies seems to have a large H^ radial variation. 
NGC 1351 A (the first dark blue arrow on the left side of 
Figure [TOl l has amongst the lowest number of bins available 
for the fit of its gradients. Despite an apparently impressive 
change in H^, careful consideration of the uncertainties forces 
us to discard any significant gradient in this bulge. The case 
of IC 5264 (magenta arrow) is diff'erent. Its spectra have high 
signal-to-noise ratios and its radial sampling is good. As for 
IC 5176, seen previously, we indeed face a case of younger 
ages in the outer parts of the bulge. But again, we do not find 
any straightforward explanation for this galaxy to be different 
from the bulk of our sample. 

The three remaining galaxies, NGC 522, NGC 1886, and 
ESO 443-042 have null or very weak gradients. They are all 
well sampled in radius, therefore the absence of spectral ra- 
dial variation can be considered real and attributable to their 
small central velocity dispersion and low global chemical en- 
richment. 

6.3. Comparison with elliptical galaxies 

The work of Sanchez-Blazquez et al. (2006) gives us the op- 
portunity to compare a set of indices in elliptical galaxies to our 
data, in a range of velocity dispersion very similar to the one 
of our sample of bulges. In particular, for our present interest, 
it allows a look at the relations between the spatial distribution 
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Fig. 11. Comparison of gradients of well-measured indices for 
bulges in this paper with those for elliptical galaxies in the 
sample of Sanchez-Blazquez et al. (2006; labeled "S-B+06"). 
Error bars are only plotted if they exceed 0.03 dex in a given in- 
dex gradient. For ease of comparison, boxes with dashed lines 
are drawn in a fixed position for each pair of gradient-gradient 
plots. 



of different indices: they offer clues on the star formation histo- 
ries in the galactic systems. Figure[TT]displays this comparison 
for gradients in Mg2, C4668', Fe4383' and Hy'^. Obviously, 
the gradients in bulges agree very well with those in elliptical 



\ o - 

8 



<U O 



< 
O 





1.8 



2.2 



2.4 



Log(7o 



Fig. 12. The relation between ALog(Age), A[Fe/H] , and 
A[c!'/Fe] with the bulge central velocity dispersion. As are mea- 
sured from the bulge effective radii to the bulge centres 



galaxies, both in amplitude and in the relations between sets of 
gradients. This strengthens further the similarities between the 
two types of spheroids. 



6.4. Quantification 

We now try to quantify the [a/Fe], age and [Fe/H] variations 
between the bulge effective radii and their cen tral parts . We 
use the models of Thomas et al. (|2003[ |2004|) . First, we de- 
rive [a/Fe] from the Mgb'-<Fe>' diagram. Subsequently and 
independently, we calculate the age and metallicity, as infeiTed 
from the H^-Fe4' plane, at r^yy and in the bulge central parts. 
The ages and metallicities are derived at the computed [a/Fe] 
values. For those, we have considered initial mean ages of 10 
and 5 Gyr, and checked that the dependence of the final age and 
metallicity variation on this a priori choice was negligible. 

The results are listed in Table H] As noticed previously in a 
qualitative way, the outer parts of the bulges have (with a few 
exceptions) higher [a/Fe], older ages and lower metallicities 
than the inner parts. Only two galaxies, IC 5264 and IC 5176, 
have significant positive age gradients. Note that due to the un- 
even spacing between the model lines in the H^-Fe4' diagram, 
the absolute errors attached to variations in age are larger when 
measured for the oldest populations (for which model lines 
are close together) than for younger populations. Similarly, the 
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Fig. 13. Generalized histograms of the distribution of radial 
changes in [Fe/H] and Log(age), calculated between the bulge 
effective radii and the bulge centre, as tabulated in Table|4] The 
dotted lines correspond to the full sample of galaxies; the solid 
lines show the histograms of the galaxies with strong gradients. 

variations themselves are maximized. Therefore, we are effec- 
tively measuring upper limits of Aage. This effect does not ex- 
ist for [Fe/H] for which the grid lines are more evenly sepa- 
rated from one another The uncertainties in A[Fe/H] are mainly 
linked to the errors in A[a/Fe] (besides the observational ones). 
This is particularly true for galaxies which seem to harbor a 
positive metallicity gradient. Galaxies having insignificant age 
gradients also have negligible metallicity gradients. Figure [T2l 
illustrates the results of Table |4] and presents the relations be- 
tween the variations in [a/Fe], age and [Fe/H] with the bulge 
central velocity dispersion. The description of the Mg2 - (To 
relation in Sect. 5.2 is applicable here as well. This is partic- 
ularly true for A[Fe/H] which has the largest dispersion: I.e., 
there isn't any significant correlation between metallicity gradi- 
ent and velocity dispersion, strictly speaking. Instead one sees 
a decrease in the range of possible gradient values at lower cr. 

In order to quantify the qualitative statement made earlier 
that bulges have comparable properties at their effective radii, 
while they differ in their central regions, we measured their 
mean ages and metalhcities. At identical ages, the weak gra- 
dient bulges have an error weighted mean [Fe/H](reflF)= 0.11 
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-0.38+ 0.31 


ESO079-003 


0.0+ 0.2 


-1.1+6.3 


-0.12+ 0.35 


ESO234-053 


0.1+0.1 


3.0+ 3.7 


-0.12+ 0.43 


ESO311-012 


-0.1+0.1 


-1.1+ 1.5 


-0.08+ 0.12 


ESO443-042 


0.0+ 0.1 


0.7+ 0.9 


-0.42+ 0.27 


ESO512-012 


0.1+0.2 


-0.1+7.2 


0.00+ 0.77 



Table 4. Variation in [a/Fe], age and [Fe/H] between the bulge 
bulge effective radii and the central regions. Negative values 
indicate an increase of the considered quantity. 



dex, with a dispersion of 0.14 dex. Strong gradient bulges have 
a mean [Fe/H](reff)=0.05 dex with a dispersion of 0.19 dex. 
This must be compared to a mean central metallicity of [Fe/H] 
= 0.13 dex (dispersion of 0.19 dex) and 0.45 dex, (dispersion 
of 0.14 dex), for the weak and strong gradient bulges, respec- 
tively. Note that these values are derived from SSP models. As 
such, they are not definite. The comparison they allow between 
galaxies is however robust. The galactocentric distance of one 
effective radius is chosen because all our galaxies have spectra 
of adequate quality there. However, it does not yet sample the 
most outer bulge regions, where we think there is a lot to learn 
on the very early stages of bulge formation. 

Figure [13] shows the generalized histograms of ALog(age) 
and A[Fe/H]. We choose here a logarithmic representation of 
the age for internal homogeneity (in the sense that index vari- 
ations scale approximately linearly with Log (age)) and for a 
more direct comparison with the metallicity scale. Both the 
full sample, in dashed line, and the strong-gradient galaxy sub- 
sample, in solid line, are shown. By nature of the use of gener- 
alized histograms, the errors are taken into account. 

Both the age and the metallicity gradient distributions are 
very well peaked, at ~0.15 dex in Log(age) (1.5 Gyr in age) 
with a dispersion cr = 0.2 dex (in Log(age)) (1.3 Gyr in age), 
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and at -0.4 dex with a dispersion cr - 0.3 dex for [Fe/H]. They 
correspond to peak values of the gradients of ~0.07 and ~ -0.2 
for the Log(age) and [Fe/H], respectively. 

Similarly, the mean gradients for the age and [Fe/H] are 
0.06 and -0.16, respectively. In other words, metallicity gra- 
dients are two to three times larger than those in age (in log 
scale). These values are very close to those derived for ellip- 
tical galaxies with a comparable range of velocity dispersions 
by Sanchez-Blazquez et al. (12006 , see summary of other anal- 
yses therein): 0.082 and -0.206. The different grids of mod- 
els and relative distribution in velocity dispersion within each 
sample suffice in explaining the small differences between their 
and our analyses. Although the distribution of gradient ampli- 
tude is more meaningful than their mean values in studying the 
processes at play in building spheroids, the similarity between 
mean values of gradients in ellipticals and bulges suggests that 
they must share a fair amount of common formation history. 

7. Conclusions 

We have presented the analysis of radial gradients of stellar ab- 
sorption line strengths in a sample of 32 edge-on spiral galaxy 
bulges. The sample galaxies span nearly the full Hubble se- 
quence (from SO to Sc types), and have a large range of dy- 
namical properties, with their bulge central velocity dispersions 
ranging from ~ 60 to 300 km s ' . Our main conclusions are the 
following: 

• Most bulges do present radial stellar population gradients. 
The outer parts of bulges harbor weaker metallic absorption 
lines than the inner regions. The distribution of these gradients 
among bulges are generally well peaked. They also display a 
real intrinsic dispersion, implying the presence of a variety of 
star formation histories within a common framework. 

• In a number of cases, the gradients in Mgi and Mg2 do not 
follow those in Mgb. We suggest that this is due to the inclu- 
sion of carbon, traced by the C4668' index, in their bands. This 
explains a number of apparent previous discrepancies between 
works dealing either with Mg2 or Mgb, in particular for ellipti- 
cal galaxies. 

• We argue that the existence of a correlation between gra- 
dients and central velocity dispersion among bulges is not an 
appropriate terminology. Instead, one sees that bulges with 
large velocity dispersion can exhibit strong gradients (they can 
also have negligible ones), while this probability diminishes at 
lower cr. Below 125 km s"', we do not find any bulge with sig- 
nificant radial spectral gradient. The same kind of dual regime 
with the velocity dispersion was observed in elliptical galaxies. 
This gradual build-up of the index vs. cr relation can only be 
clearly observed for indices with large dynamical ranges. 

• The Hubble type of the parent galaxy appears to be a sec- 
ondary parameter, earlier types having statistically larger ve- 
locity dispersions than later ones. The strength of the gradients 
depends only weakly on the bulge effective radius. The differ- 
ence in sensitivity to the effective radius on one hand, and to the 
velocity dispersion on the other hand, suggests that the depth 
of the gravitational potential in which bulges are embedded is 
the main property tracing the spatial distribution of their stellar 
population. 



• Strong gradients are found in bulges with the most metal-rich 
central regions, while bulges with strong and weak gradients 
have comparable properties at the bulge effective radius. This 
indicates that the external regions of the bulges, which are typ- 
ically the oldest and most metal-poor and therefore the first to 
form, share some universal properties. This also argues in fa- 
vor of an outside-in scenario for the formation and evolution of 
bulges. 

• The analysis using SSP models indicates that radial variations 
in luminosity-weighted mean metallicity are twice to three time 
as large (in logarithmic scale) as the variations in age. While 
[Fe/H] at the bulge effective radii are on average 0.4 dex (cr = 
0.3 dex) lower than in the bulge central regions, the age dif- 
ference is of the order of 1.5 Gyr (cr = 1.3 Gyr), the inner re- 
gions being younger We have found only two galaxies with 
convincing signs of an 'inverted' age gradient. The changes 
in [a/Fe] are small (of the order 0. 1 dex) and rather constant 
among bulges. 

As to the fundamental question of the influence of the disk 
on the bulge evolution, we make the following comments: In 
terms of stellar population properties, we do not find any ob- 
vious differences between ellipticals and bulges. On the con- 
trary, it seems that our observations strengthen their resem- 
blance. One requirement to properly compare the two families 
of spheroids is to match systems of similar velocity dispersion. 
The imprint of the disk influence is generally seen via t he pres- 



ence of a bar in spiral galaxies. ^ Bournaud et al. ( l2005i) deter- 



mined that a cycle of bar formation and dissolution takes about 
2 Gyr for Sb-Sc galaxies with masses and radii comparable to 
that o f the Milky Way. Gas accre tion allows the bar to reap- 
pear d Bournaud & Combes 2002). In terms of time scales, a 
scenario where bulges would be formed through secular evolu- 
tion could thus be conceivable for late-type galaxies. However, 
for the time being, our study introduces two severe constraints 
to this type of scenario: (i) The action of the bar must preserve 
the observed gradients and star formation time scales (as mea- 
sured by [a/Fe], whose measured values are significantly larger 
than those of disk stars); (ii) The age gradients produced by the 
bar cycles (formation to dissolution) must be consistent with 
those measured (peaked at ~ 1 .5 Gyr). For example, if 2 to 3 bar 
cycles are necessary to build the bulge of an early-type galaxy 
from disk material, the bulge assembly timescale is larger than 
the one implied by the observed gradients. 

We note that these remarks also hold for models where 
bulges grow by accretion of sa tellites, combined with th e 
accretion of disk material (e.g., lEliche-Moral et al. 



Nevertheless, bars are ubiquitous among spirals, and their pres- 
ence and action could explain part of the dispersion seen among 
the bulge properties at fixed velocity disp ersion, in particula r 
by wiping out strong gradients (see, e.g. iFriedli et alJ 1994 ). 
Contrary to what was found by Moorthy & HoltzmanI ( 2006 ), 
we do not find that younger ages in the centres of bulges are re- 
stricted to barred and/or boxy/peanuts shape galaxies. In fact, 
we do not find any difference between the stellar properties of 
barred and non-barred galaxies in our sample, as seen along 
the ir bulge minor axi s. 



Kobavashil ( 120041) has proposed, as a result of her simula- 



tions of eUiptical galaxies, that galaxies with strong gradients 
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originate from small and rapid mergers mimicking an initial 
(monolithic) collapse, whereas galaxies with weak gradients 
would reveal a sequence of more significant mergers (i.e., with 
small mass ratios). In a scenario where bulges and ellipticals 
form in a similar way, i.e., disks settling after the bulk of the 
stellar population in bulges is in place, this hypothesis would 
also fit part of the observations presented here. However, bulges 
with low velocity dispersions would then require a different for- 
mation mechanism. 

Acknowledgements. This work was supported by the Spanish research 
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